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ABSTRACT 

Two scenarios of possible ion heating due to finite amplitude parallel propagating Alfven waves in 
the solar atmosphere are investigated using a ID test particle approach. 1. A finite amplitude Alfven 
wave is instantly introduced into a plasma (or equivalently, new ions are instantly created). 2. New 
ions are constantly created. In both scenarios, ions will be picked up by the Alfven wave. In case 
1, the wave scatters ions in the transverse direction leading to a randomization (or heating) process. 
This process is complete when a phase shift of ±7r in the ion gyrospeed is produced between particles 
with characteristic parallel thermal speed and particles with zero parallel speed. This corresponds to 
t = -j^- {k is the wavenumber and v t h is the ion thermal speed) . A ring velocity distribution can be 
produced for a large wave amplitude. The process yields a mass-proportional heating in the transverse 
direction, a temperature anisotropy and a bulk flow along the background magnetic field. In case 2, 
continuous ion creation represents a continuing phase shift in the ion gyrospeed leading to heating. 
New particles are picked up by the Alfven wave within one ion gyroperiod. It is speculated that the 
mechanism may operate in the chromosphere and active regions where transient events may generate 
finite amplitude Alfven waves. To appear in ApJ Letters May 10, 2007. 
Subject headings: Waves; Sun: chromosphere; Sun: corona 



1. INTRODUCTION 

It has been widely speculated that the energy that 
heats the corona comes from the convective flows in the 
photosphere. The energy is somehow transported into 
the coronal part through the magnetic field. It is natu- 
ral to think that Alfven waves channel the energy to the 
corona. Indeed, these waves have been observed in the 
solar atmosphere (Ulrich 1996), and are ubiquitous in 
the extended corona — the fast solar wind (Smith et al. 
1995). However Alfven waves are difficult to dissipate in 
collisionless plasmas. Hence, non-linear processes have 
been assumed to cascade the wave energy from low to 
high frequencies where wave dissipation is readily possi- 
ble (Hollweg 1986, Li and Habbal 2003). 

When Alfven waves propagate in a partially ionized 
plasma, neutral-ion collisions produce a channel for the 
wave dissipation (De Pontieu & Haerendel 1998; De Pon- 
tieu et al. 2001; Leake et al. 2005). Ions and neutrals 
can be collisionally coupled. A slippage between the ion 
and neutral populations leads to the wave dissipation. If 
we consider neutrals are constantly ionized at the chro- 
mosphere, we will show that these newly created He +1 
ions will be picked-up by the wave and will be energized. 

In this Letter, a new scenario, ion pickup by an Alfven 
wave, is explored. The pickup process can lead to the 
heating of ions, and it must also dissipate the wave. In 
section 2, we discuss the ion pickup process by an in- 
stantly introduced wave (or identically, the ions are cre- 
ated instantly). In section 3, the pickup of continuously 
created ions by an Alfven wave is investigated. Finally, 
in section 4 we discuss possible applications of the ion 
pickup process in the solar atmosphere. 

2. ION PICKUP BY ALFEN WAVES 

Consider a parallel propagating monochromatic disper- 
sionless Alfven wave with angular frequency to and wave- 



number k, and to — kvA(vA is the Alfven speed), 
wave electromagnetic field 5B W and iE M are 



5B W = Bk{cos<f>kix - sin< 
5~E W = -v A B a /B x 6B n 
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and iy are unit vectors, 



Bq the background magnetic 
field, and <pk = k{v A t — z) denotes the wave phase. The 
motion of a particle is described by 
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where £Iq = ejBo/rrij, £l k — £jBk/rrij, and j refers phys- 
ical quantities of ion species j. As a first order approx- 
imation, t>|i w W||(0) is a constant, where vn(0) is the 
initial ion parallel velocity. The approximation is valid 
when r2fe / O = B^/Bq is small and the wave frequency 
is low so |f?o| 3> |A;(ui| — va)\. With the initial condition 
uj_ = uj_(0) and z — z(0), the solution of Eq. ^ for a 
low beta plasma is (Wu et al. 1997, Lu et al. 2006) 
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where z — z(0) + v^(0)t and we have adopted approx- 
imations flo — k[vA — W||(0)] W fi(); V A — W||(0) w va- 

The last term of u± is the perturbed ion velocity of an 
Alfven wave. The first term is due to the gyromotion of 
the particle and the modification of the gyromotion due 
to the wave. Note, if uj_(0) = -v A {B k /B )e- ikz( -°'> and 
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ui|(0) = 0, the particle is surfing on the wave, there is no 
gyromotion. 

Now let us consider an ensemble of newly created par- 
ticles with a Maxwellian distribution. The particles' av- 
erage parallel speed is zero but the transverse component 
is u±f 

Bi. , 

«±(0) = u± r (0) + u_l/(0) = u± r (0) - av A -fe ikz , (6) 

-DO 

where a, a constant, describes the degree that the newly 
created particles are settled in the wave field. If a = 
1, the newly created particles are already settled. uj_ r 
denotes the random perpendicular thermal speed. After 
time t particles with different initial position z(0) = z — 
vn(0)t will arrive at z. The average transverse velocity 
at z is 



B k 

U± = -v A —e 
Bo 
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where A k = -4= f°° cos(kvthtx)e x dx = e 

1/2 



I'll, 



( 2k B Tj \ 
V m 3 ) 



and Tjo is the initial temperature 

of species j. Eq ((7|) illustrates the pickup of ions by the 
wave. Note Ak « when kv t ht > n. Hence, regardless 
of the value of a, these particles will be picked up by the 
Alfven wave: eventually U± is determined by the wave 
only. Subtracting |(7J) from ([5]), one finds the random 
perpendicular velocity at z: 

u±-U ± = u ±r (0)e- mot + (1 - a )v A —e lk[z - v « <- ^ e ~ inot 

Bo 



(1 - a)A k v A 



-°fc A(kz-Uot) 
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Eq. (8) describes the particles' gyrospeed in the frame 
of U±. When a = (initially the newly created par- 
ticles have zero flow speed), particles will be strongly 
scattered in the phase space. The second term on the 
right hand side is crucial for our understanding the heat- 
ing in the pickup process. In the phase space frame 
of U±, a randomization process is complete when parti- 
cles with characteristic speed v» (0) = ivth moving from 
z(0) = z — f||(0)i to z are scattered in the perpendicular 
direction and a phase shift ±7r in the gyrospeed relative 
to particles with v\\ (0) = is produced. This translates 
to kvtht ~ 7T. For a finite amplitude Alfven wave, if 
v A Bk/Bo > Vth, the distribution function will be a ring. 
Hence particles are strongly heated. On the other hand, 
if a = 1, initially particles are already "picked up" by 
(surfing on) the wave. Only the initial random motion 
survives, and there is no heating. One notices that if 
vth = 0, Ak = 1, there is also no heating. Hence, the 
heating is a warm plasma effect. With U±, the perpen- 
dicular temperature can be found 
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where f3j = 8im e kBTjo/BQ is the plasma beta of species 
j. The heating is obviously mass-proportional. 

Substituting Eq. (O into Eq. (0}, and letting v A — 
f ii (0) ss v A and Slo — k(v A — i>ii(0)] ~ f2o, one finds 

B 2 

v\\ = U||(0)+i; A (l-Q;)^|{l-cos[fioi-A;^t-feu||(0)t]}, 
a o 

(10) 

Using the same procedure to obtain U± and T±, we find 
the average parallel velocity and temperature: 

B 2 

u \\ =v A -£(l-a)[l-A k cos(n t-kv A t)}, (11) 
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m p /3jB 2 

cos 2 (fV - kv A t) + (1 - C fe ) sin 2 (fV - kv A t)}] , (12) 

where C k = ;4= J^ ^ cos 2 (kvthtx)e~ x dx = 0.5 + 

0.5e- fc2t 't 2 '>* 2 . When i -> oo, A fc = 0, and C k = 0.5. 
We then obtain asymptotic values: 

B 2 

So 
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Subtracting (fTTj) from (|10[) . when £ — > oo, one finds 

i? 2 

W|| — U\\ = v^(0) + v A (l — a)—gcos[Qot — kv A t — kv\\(Q)t\. 

''a 

(15) 
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Fig. 1. — The time history of (a) the perpendicular temper- 
ature T±ff e /TneO and t ne parallel temperature T^ H< ,/T[] e Q, (b) 
the temperature anisotropy Tj_ He /T^ He , and (c) the average par- 
allel velocity Un/v A . Here we choose the particle as He +1 , and the 
parameters are: a = 0, u/(l p = 0.25?r/125 » 0.0063, /3 He = 0.01, 
B V B = °- 09 - 
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Fig. 2. — Top panel: Scatter plots of helium He +1 velocity in the velocity plane perpendicular to Vi\. Bottom panel: scatter plots of 
helium He^ 1 velocity in the plane parallel to the background magnetic field. Relevant parameters are those in Fig.l. 

the frame of the Alfven phase speed, the kinetic energy 
of particles is conserved. As particles are picked up by 
the wave, they gain a bulk acceleration along Bo- 

3. PICK-UP OF CONTINUOUSLY CREATED IONS 

In a partially ionized plasma, the creation of new ions 
due to photoionization or collisional ionization is highly 
random. Consider a region with a low plasma beta and 
where the amplitude of an Alfven wave is finite. When 
ions are created from neutrals, they are subject to rapid 
gyromotions. The continuous creation of ions represents 
a continuous phase shift in the ion gyrospeed. This will 
naturally mix the phase of particles' gyrospeeds. To ex- 
plore the pickup process of continuously created ions, a 
new test particle simulation is conducted. The Alfven 
wave and the simulation box are the same as in section 
2. However, 1000 particles are deployed every 0.2517^ 1 
at < flpt < 100 (one particle per cell). The beta value 
of the newly created He +1 ions is 0.01 and they have 



It is possible that the second term on the right hand 
side is larger than the initial random thermal velocity 
f|| (0). In this case, the velocity distribution plotted in 
the up — v y plane may be a ring as well (see Fig. 2). 

To verify the above analysis, a test particle simu- 
lation was conducted by calculating the full dynamics 
of particles (Eq. \2§ in the electromagnetic field of a 
given Alfven wave described by Eq. 1. The equations 
were solved using the Boris algorithm with time step 
At = 0.01f2~ 1 (f2 p is the proton gyrofrequency) and peri- 
odic boundary conditions. Initially, 200000 particles with 
Maxwellian velocity distribution were evenly distributed 
in a region with length 1000^0- 1 (1000 cells). The 
flow speed of these particles was zero. The average par- 
allel velocity, and the parallel and perpendicular temper- 
atures were obtained by using the following procedure: 
we first computed U\\ = (v z ), T\\ = j£ ((v z - (v z )) 2 ), 

and T± = ^((v x - (v x )) 2 + (v y - (v y )) 2 ) in every cell 

(() denotes an average over a cell). These quantities were 
then averaged over all cells, so only random motions con- 
tribute to the temperatures. Fig.l shows the time evo- 
lution of (a) the perpendicular temperature T\\ He /TH e0 
(here refers to the initial value) and the parallel tem- 
perature T±He/Tne0, (b) the temperature anisotropy 
Tj_/Tj|, and (c) the average parallel velocity Uii/va of 
the test particle simulation. The "shaded areas" are due 
to rapid ion gyrations. The results are consistent with 
our analytical predictions. At £l p t f=a 12000, an asymp- 
totic stage is reached: Uu/va ~ 0.086, T±jj e /Tjj e Q w 34 
and Ti\H e /TH e o w 3.9 with the temperature anisotropy 
Tj_He/T\\ He ps 8.7. 

Fig. 2 displays scatter plots of the helium ions in the 
five cells around z = 500v A ^p 1 at n p t = 0, 2800, 6600 
and 11000. Initially, helium ions satisfy the Maxwellian 
distribution with vth = 0.005UA and fine = 0.01. Helium 
ions are dramatically scattered in the transverse direction 
in the phase space. The ring velocity distributions are 
nicely shown in Fig. 2. The ring velocity distribution is 
unstable (Lu and Wang 2005). The velocity distribution 
may eventually be thermalized by relevant microscopic 
instabilities of the ring distribution. The bulk accelera- 
tion of particles along Bq is obvious. This is because in 
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Fig. 3.— The time history of (a) T ±He /T He0 and T^ He /T He0 , 
(b) Tj_ He /T|| He , and (c) U\\/v A . Here kv A /£lp = 0.25tt/125 w 
0.0063, f3 He = 0.01, Bl/Bl = 0.09. The He +1 ions are steadily 
and uniformly created between 2 = and 1000^0^ at < f2 p t < 
100. 
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Fig. 4. — Top panel: Scatter plots of helium He + velocity in the velocity plane perpendicular to Vu . Bottom panel: scatter plots of 
helium He^ 1 velocity in the plane parallel to the background magnetic field. Relevant parameters are those in Fig. 3. 



zero flow speed. The results of the simulation are shown 
in Figs. 3 and 4, which are plotted in the same way as 
Figs. 1 and 2. Note that in the v x — v y plane, because 
the radius of the ring (Fig. 3) is BkVA/B , the origin, 
where new particles are created around there, is always 
on the ring. When new He +1 ions are continuously cre- 
ated, the heating and scattering process are completed 
when they have filled the ring of radius B^va/Bq in the 
v x — v y plane. This only takes one gyroperiod of He +1 
ions, or when Vt p t = 87r, as shown in Fig. 3 and 4. Hence 
the heating of these particles is extremely rapid. The 
ion creation process was terminated at time t = lOOfi" 1 , 
however no change is observed as Vl p t > 87r. We believe 
that the rapid heating due to Alfven waves reported by 
Wang et al. (2006) is due to their deployment of parti- 
cles within one gyroperiod, even though they used many 
wave modes in their computations. 

4. DISCUSSION FOR POTENTIAL APPLICATIONS 

The Alfven wave pickup of newly created ions may 
have applications at the top of the chromosphere or lower 
transition region. Let's consider the region with an ion 
cyclotron frequency much higher than Coulomb collision 
frequencies, low plasma beta and finite amplitude Alfven 
waves. For neutral helium, when wave frequencies are 
higher than the ion neutral collision frequencies, the per- 
turbed ion and neutral velocity will have a considerable 
difference between their phases and amplitudes (De Pon- 
tieu & Haerendel 1998). In coronal funnels, the plasma 
beta at the lower transition region is much smaller than 
unity (Li 2003). For finite amplitude waves, the pickup 
process described above may operate and heat newly cre- 



ated He +1 ions. Due to collisions, the temperature is 
not expected to reach the value given by (fl4|) . Instead, 
the energy will be transferred to protons and electrons 
via Coulomb coupling. Obviously, the thermal energy of 
He +1 gained through the pickup process will be passed 
to He +2 later on when He +2 are created. As the ioniza- 
tion continues to the transition region (Hansteen et al. 
1997), the pickup process will contribute to the plasma 
heating there. 

The pickup of instantly created ions by an Alfven wave 
or the pickup of ions by an Alfven wave instantly intro- 
duced into a plasma is quite efficient even for low fre- 
quency waves in collisionless plasmas. The pickup pro- 
cess is complete when kvtht = 7r or t = ^j^tA, where 
tA is the period of the wave. Transverse oscillations of 
active region coronal loops that suddenly appear with 
periods of 2—33 minutes have been frequently observed 
by TRACE (Nakariakov et al. 1999; Aschwanden et al. 
2002). The oscillations usually just last several wave pe- 
riods. The relative amplitude B^J B of these waves can 
reach 0.05 (Aschwanden et al. 2002). For typically ob- 
served 10 6 K coronal loop with density of 10 9 cm -3 and 
magnetic field of 30 G, the plasma beta value is 0.00385. 
Due to the pickup process, the perpendicular tempera- 
ture of helium ions and oxygen ions may be increased to 
3.6 and 11.4 times of their original values, respectively. 
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